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3.  RESEARCH  ACCOMPLISHED  AND  PAPERS  PUBLISHED 


The  research  accomplished  deals  with  two  main  topics: 

(a)  stretched  atoms  in  collisions  with  other  species  (electrons,  ions  and  atoms) 
and 

(b) stretched  (excited)  atoms  resulting  from  electron  -  ion  recombination.  The  dynamics  of 

stretched  atoms  has  therefore  been  concerned  with  collisions  and  with  recombination. 

The  following  research  in  these  two  topics  has  been  written  up  and  published. 

1.  Recombination  at  Ultra-low  Energies ,  by  M.  R.  Flannery  and  D.  Vrinceanu,  in  APS 
Topical  Conference  on  Atomic  Processes  in  Plasmas,  eds.  E.  Oks  and  M.  S. 
Pindzola  (AIP  Press,  1998). 

2.  The  Classical  Atomic  Form  Factor ,  by  D.  Vrinceanu  and  M.  R.  Flannery,  Physical 
Review  Letters,  (1999) 

3.  Cross  Sections  for  Electron  Excitation  of  the  2 3S  Metastable  Level  of  He  into  Higher 
Triplet  Levels ,  G.  A.  Pech,  M.  E.  Lagus,  L.  W.  Anderson,  C.  C.  Lin  and  M.  R.  Flannery, 
Phys.  Rev.  A  55,(1997),  pp.  2842-2856. 

4.  Elastic  Scattering:  Classical,  Quantal  and  Semiclassical,  by  M.  R.  Flannery  in 
Atomic,-  Molecular  and  Optical  Physics  Reference  Book,  ed.  G.  W.  F.  Drake 
(AIP  Press,  New  York,  1996),  Chap.  43,  pp.  499-525. 

5.  Electron-Ion  and  Ion-Ion  Recombination,  by  M.  R.  Flannery,  in  Atomic,  Molecular 
and  Optical  Physics  Reference  Book,  ed.  G.  W.  F.  Drake  (AIP  Press,  New  York, 
1996),  Chap.  52,  pp.  605-629. 

The  first  two  papers  are  attached  as  Appendices  A  and  B  of  the  present  report.  Publication 
#3  represented  a  joint  collaboration  with  the  experimental  collision  group  headed  by  Pro¬ 
fessor  C.  C.  Lin  at  the  University  of  Wisconsin  on  electron  -  excited  atom  collisions.  The 
first  page  of  this  publication  is  reproduced  in  Appendix  C.  Publications  #4  and  5  were  fully 
documented  in  the  final  report  AFOSR  Grant  F  49620-94-1-0379,  but  were  only  published 
during  the  current  research  period. 
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4.  INVITED  PAPERS  PRESENTED  AT  CONFERENCES 

The  following  papers  were  presented  at  professional  scientific  conferences  during  the  period: 

1.  Recombination  at  Thermal  and  Ultracold  energies ,  M.  R.  Flannery,  M.  R.  Flannery, 
Invited  Paper,  presented  at  50th  Gaseous  Electronic  Conference,  Madison, 
Wisconsin,  October  5-9,  1997. 

2.  Recombination  at  Low  and  Ultralow  Temperatures ,  M.  R.  Flannery,  M.  R.  Flannery, 
Invited  Paper,  presented  at  Molecular  Ion  Physics  Workshop,  Oak  Ridge  Na¬ 
tional  Laboratory,  Oak  Ridge,  Tennessee,  February  20-21,  1998. 

3.  Recombination  at  Ultracold  Energies,  M.  R.  Flannery,  Invited  Paper,  presented  at 
Atomic  Processes  in  Plasmas,  11th  Topical  APS  Conference,  Auburn,  Alabama, 
March  22-26,  1998. 

4.  Recombination ,  M.  R.  Flannery,  Allis  Prize  Invited  Plenary  Lecture,  presented 
at  American  Physical  Society  (APS)  Division  of  Atomic,  Molecular  and 
Optical  Physics,  Santa  Fe,  New  Mexico,  May  27-30  ,  1998. 

5.  Three-Body  Recombination,  M.  R.  Flannery,  Allis  Prize  Invited  Plenary  Lecture, 
presented  at  51st  Gaseous  Electronics  Conference,  Maui,  Hawaii,  October  19-22, 
1998. 

6.  Three-Body  Electron-Ion,  Ion-Ion  and  Neutral-Neutral  Recombination  at  Thermal  and 
Ultra-Cold  energies,  M.  R.  Flannery,  Invited  Paper,  presented  at  65th  Meeting  of 
the  Southeastern  Section  of  The  American  Physical  Society,  Miami,  Florida, 
November  13-15,  1998. 
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5.  HONORS  AND  AWARDS  RECEIVED 


The  following  academic  and  scientific  Honors  and  Awards  were  received  during  the  period: 

*  Doctor  of  Science  ( honoris  causa):  The  Queen’s  University  of  Belfast,  1998 
Cited:  For  distinction  as  a  physicist. 

*  Recipient  of  the  APS  Allis  Prize  for  the  Study  of  Ionized  Gases,  1998 

Cited:  For  advancing  the  understanding  of  recombination  processes,  in  particular 
for  developing  a  microscopic  theory  of  three-body  ionic  recombination;  and  for  his 
novel  applications  of  classical  and  quantum  mechanical  methods  to  the  dynamics  of 
atomic, molecular  and  ionic  systems. 

*  Elected  Honorary  Member:  Royal  Irish  Academy,  1997 

Cited:  In  recognition  for  his  distinguished  services  in  the  Section  of  Science. 
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Appendix  A 

Recombination  at  Ultra-low  Energies 

M.  R.  Flannery  and  D.  Vrinceanu 


School  of  Physics,  Georgia  Institute  of  Technology,  Atlanta,  GA  30332-0430 1 


Abstract.  Three-body  electron-ion  recombination  is  described  at  ultralow  electron 
temperatures  Te.  At  4K,  the  initial  stage  involves  extremely  rapid  collisional  capture 
into  high  Rydberg  states  n  >  200  with  high  angular  momentum  Iwn-lata  rate 
~  T“4  5.  This  is  followed  by  extremely  slow  collisional-radiative  decay.  The  key 
collisional  mechanism  appears  to  be  collisional  /-mixing  of  the  Rydberg  atoms  A(n)  by 
ions  and  electrons  until  sufficiently  low  Ts  are  attained  so  as  to  permit  relatively  rapid 
radiative  decay  to  the  lowest  electronic  levels.  This  sequence  is  in  direct  contrast  to  the 
sequence  of  much  slower  collisional  capture  at  higher  Te  followed  by  the  much  faster 
decay  of  A(n)  by  electron  collisions  to  lower  levels  where  radiative  decay  completes 
the  recombination.  At  ultra-low  temperatures,  the  rate  limiting  sequence  is  therefore 
collisional  /-mixing  followed  by  radiative  decay  in  contrast  to  recombination  at  much 
higher  energies  and  electron  densities  Ne  ~  108  cm-3,  where  the  rate  limiting  step  is 
the  initial  collisional  or  radiative  capture  at  intermediate  Te(~  1  eV)  and  higher  Te 
(~  10  eV)  respectively. 

An  exact  classical  solution  for  /-mixing  is  obtained. 


INTRODUCTION 

Three-body  collisional  capture 

e“  +  Az+  +  e-  —  A^+(n)  +  e"  (1) 

and  radiative  capture 

e-  +  Az+— »A  (z~V+(n)  +  hu  (2) 

in  an  electron  gas  at  cryogenic  temperatures  T  ~  4 K  and  lower  are  considered. 
Simple  classical  expressions  for  the  cross  section  a  and  rate  a  ~  Z3T~4  5  of  (1) 
are  derived.  We  shall  also  provide  an  exact  classical  solution  for  the  angular  mo¬ 
mentum  mixing  n,l  — ►  n,  V  of  Rydberg  levels  by  the  time  varying  electric  field 
£  =  -[Zie/R3(t)]R(t)  generated  by  distant  collisions  with  an  ion  of  charge  Z\e. 
The  work  has  particular  significance  to  the  production  of  anti-hydrogen  at  4  K  and 
to  proposed  experiments  on  electron  recombination  with  Rb+  and  Xe+  at  5  mK. 


*)  This  Research  is  supported  by  the  U.S.  Air  Force  Grant  No.  F49620-96-1-0142 


RADIATIVE  RECOMBINATION 


The  quantal  cross  section  for  photoionization  of  A(n ,  /),  in  conventional  notation 
[1],  is 

/ 1  \  f  64o  .  ( In  \31  21  +  1  ^  .  .  .  .2/4-1  .  . 

<T'('“/)  =  U3V5(™o)z]  (w  }  — G«lhv) 

where  k&n(his)  is  the  semiclassical  (Kramers)  photoionization  cross  section  [2]  and 
the  bound- free  Gaunt  factor  Gni  ~  1  and  as  u  =  hv/In  tends  to  1  and  oo 

respectively.  The  former  limit  is  appropriate  to  ultracold  ejected  electrons  while 
the  latter  limit  indicates  that  low  /'  are  favored.  The  detailed  balance  relation  is 

9e9A^R(E)  =  gvgAkl(Ji{hu) 

where  the  statistical  weights  are  (2, 1, 2, 2(2 1  +  1)}  for  e,  A+,  hu  and  A,  respec¬ 
tively  and  where  k%/k2  =  a2(hu)2/(4EIH)  =  (v0jc)2(IH/4E)  with  Ih  =  (e2/2ao). 
The  cross  section  for  radiative  recombination  of  electrons  with  energy  E  is  then 

,  x  \{Sna2a3\  2Z2I„\  [2  In  1  „  /r,  ,  r  N 

^  =  U-S^rJ  —  I  ln(vTB)J  +  « 

which  is  much  smaller  by  a  factor  0(a2)  than  <7/.  Gn  is  the  /-averaged  Gaunt  factor 
n~2  E"=o1(2/  +  1  )Gni.  The  corresponding  rate  is 


&RR(n,  Te )  =  ( ve ) 


8™2a3\  [2  Z2IH] 

-jtt  -w  nFM 


-  „  j. 


Fjb„) = b,r  J 


where  bn  =  In/kTe  and  u  =  hv/In.  For  ultracold  recombination,  Gn  ~  1  so  that 


<^RR{n,  Te  — *  0) 


-M'S 


2 &1h 
kTe 


-bnebnEi(bn) 


=  OoiTJ-h^ibn) 

Tl 


the  semiclassical  (Kramers)  rate  [1]  where 


o0(Te)  =  5.966  1(T13  (^)V2Z 


2  cm3s-1 


and  where  Ei  is  the  first  exponential  integral  /b°°  e~x  dx/x.  Since  bebEi  increases 
from  0  to  1  as  b  increases  from  0  to  oo,  (3)  shows  that  Radiative  Recombination 
favors  transitions  to  low  n. 


Ultracold  Radiative  Recombination  into  the  high  level  n0  determined  by  kTe  =  In 
then  proceeds  at  the  rate 

aim(n )  =  7.12  10  (^— J  —  cm3s 

The  rate  into  all  levels  n'  >  n  is 

6n 

OiRR{n'  >n)  =  a0{Te)  J  ebEi(b)  db 

o 

=  Oio(Te)[^  +  ln6„  +  e6nEi(6n)] 
where  7  =  0.5772  is  the  Euler’s  constant.  For  bn  =  1,  then 

otRR{n'  >  no)  =  7  10  13  —'j  cm3s  1 

The  rate  into  all  levels  n'  <n  including  the  effect  of  cascades  [3]  is 

/  QflflX  0  63 

ctRR{n'  <  n)  =  3.80  10~12  J  Z2  cm3s-1 


COLLISIONAL  CAPTURE  RATE 


The  forward  recombination  and  reverse  ionization  rates  aR  and  a/  of  process  (1) 
are  interconnected  by  detailed  balance 


NeN(n)ai(n)  =  NeN+aR(n ) 


(4) 


where  the  Saha-Boltzmann  distribution  N  of  bound  levels  n  in  collisional  equilib¬ 
rium  with  the  continuum  state  distributions  Ne  and  N+  is 


N(n)  _  A3 
N+N0  ~  (2tt mefcTe)3/2 


n2  exp  (In/kTe)  ;  In  =  Z 


2  [h 
n2 


(5) 


where  In  is  the  ionization  energy  level  n.  It  is  useful  to  rewrite  (5)  in  the  equivalent 
forms: 


Jttt  =  z  (7tR2)3/2  (n2/n£)exp(no/n2)  ;  In/kTe  =  nl/n2  (6) 

iv+iVe 

in  terms  of  the  characteristic  Coulomb  Radius  ZRe  =  Z(e2/kTe )  =  2no(ao/Z),  and 

=  f  (nRlf/2  6n5/2  -T-r  exp  {bn)  ;  bn  =  In/kTe 
jV+ AT  2  dn 


(7) 


where  dbn/dn  =  (2 bn)/n.  The  distribution  (6)  exhibits  a  minimum  at  n  =  n0  which 
corresponds  to  bn  =  1. 

The  rate  for  ionization  of  A(n)  by  electrons  of  energy  E  is 


OO 

Qj(n)  =  (yeOj(E))  =  ( ve )  j  e  e~€  07(e)  de  ;  e  =  E/kT 

o 

A.  Detailed  Balance  Rate:  The  Thomson  collisional  cross  section  [4] 

for  electron  ejection  from  the  Rydberg  A(n),  yields  the  rate  of  ionization 
a,(n)  =  Z2tRK v,)~  (l  -  ii„e‘"E1(6n)) 

From  detailed  balance  (4)  and  (7),  the  rate  of  collisional  capture  into  level  n  is 
then: 


acc(n)  =  Z3(TRlf'2(ve)  Nt  (nd^2)"1!!  -  &„el"E,(M]  (9) 

The  factor  multiplying  [  ]  in  (9)  varies  as  Z3n42^-2.  Since  behEi(b)  — ■  1  —  ft-1  for 
very  large  b  then  acc  ~  Z3n6T’(T1  at  ultracold  energies  when  kTe  <C  In- 
Recombination  into  all  levels  E  <  n  proceeds  at  a  rate 

c*cc(£)  =  it,  acc{n')  «  j  acc{n')dri  =  jf  acc{n dbn  (10) 


where  b  =  I\/kTe  can  be  replaced  by  infinity  since  61  »  bn.  Then 


Ctec(E)  =  y  (irfi2)5/2  (»«> 


where 


Since 


6„e6nE!(6n) 

67/2 


db 


Rl/2(ve)  = 


8.094  10~9 
T(K)4b 

5.144  10-27 

EW 


=  5.769  10 


-20  (  300  \ 

\T(K)J 


4.5 


4.070  10 

EM 


-32 


for  Te  expressed  in  different  units,  then 


(11) 


(12) 


/WK4-5 

acc(£)  =  6.71  10-20JVe  (jr)  R(b „) 

which  exhibits  the  characteristic  Z3T~ 4  5  dependence  [5].  The  function  is  R(bn)  = 
I(bn)/I(  1).  Since  the  equilibrium  distribution  (7)  exhibits  a  maximum  at  kTe  = 
where  n  =  n0  or  bn  =  1,  there  is  a  bottleneck  at  n  =  no  and  it  may  be  argued 
that  all  levels  n  >  no  are  collisionally  connected  with  the  continuum  and  may  be 
regarded  as  subject  to  re-ionization  alone,  while  levels  n  <no  are  connected  mostly 
with  lower  bound  levels  and  are  therefore  subject  only  to  further  de-excitations. 
Then  R  in  (12)  is  unity. 

In  terms  of  this  bottleneck  level  n0, 

/2\ 1/2 

ZRe  -  2nl(a0/Z)  ;  (ve)  =  ( Zv0/n0 ) 

Then  (9)  provides  the  CC-rate  into  level  n0  as 

aCC(n0)  =  64  tt2  n«(a2u0)(iVeag)  {[1  -  eE^l)]  =  0.4} 


and  (11)  the  CC-rate  into  all  levels  n  <  no  as 

t *cc(n  <  n0)  =  32  7r2no(aoV0)(-^eOo)  {1(1)  =  0.1329}  (13) 

The  strong  n0- dependencies  are  therefore  highlighted.  Since  n2  ~  Z2/kTe,  the 
Z3T~i  b  dependence  of  (11)  is  recovered  from  (13).  Although  the  Thomson  cross 
section  (8)  can  be  replaced  by  a  variety  of  cross  sections  derived  from  more  elaborate 
binary-encounter  classical  theory,  the  key  dependence  Z3!1”4'5  is  still  preserved. 
Allowance  for  the  discrete  n'  summation  rather  than  the  continuous  integration  in 
(10)  shows  that  a«(n)/2  is  simply  added  to  (11). 

B.  Direct  CC  Rate:In  the  three-body  CC  (1),  some  of  the  energy  gained  by 
the  incident  (test)  electron  1  moving  with  total  energy  E\  =  (3/2 )kTe  and  kinetic 
energy  T\(R)  =  3/2kTe+\ViS(R)\  in  the  field  of  the  ion  3  with  charge  Ze  is  lost  upon 
repulsive  Coulomb  collision  with  an  ambient  (field)  electron  2  of  average  energy 
E2  =  (3/2 )kTe  in  an  electron  gas  bath  at  temperature  Te.  For  (1-2)  scattering  by 
angle  in  the  (1-2)  CM,  the  energy  in  the  laboratory  (fixed  ion  3)  frame  lost  by  1 
to  2  is 

A E  =  [Ti (R)  -  E2)  sin2  -  =  \Vl3(R)\ sin2  - 

In  order  that  the  total  energy  of  1  is  changed  from  E\  to  all  bound  energies 
E{  <  —€q kTe,  the  (1-2)  scattering  must  satisfy 


|V13(/?)|  sin2  ^  +  e0)kTe 


(14) 


In  a  repulsive  Coulomb  (1-2)  encounter 

.  2  *  _  b\ 
sm  2  62  +  bl 


bc  2  Erel 


where  Erei  is  the  energy  of  (1-2)  relative  motion  which,  on  average,  is  3/2 kTe. 
Condition  (14)  then  provides  the  averaged  cross  section 

=  '  ^  =  ^2-  Re  =  W,  (15) 

for  the  formation  of  bound  states  with  energy  E{  <  —eokT3  by  (e  -  e)  collisions  at 
test-electron  -  ion  separation  R  <  Ro-  For  R  >  -Ro,  0\2  vanishes.  When  eo  =  0,  Rq 
is  the  Thomson  trapping  radius  [1].  The  collisional  capture  cross  section  at  total 
energy  2 

E  =  \mev2{R)  +  V13{R)  =  l-meU 2  +  V13{R)  +  2  meE^ 
is 

oo  oo 

acc{E)  =  2tt  j  pdp  I  W(R,  E)  dt  (16) 

0  — oo 

where  W(R,  E)  dt  is  the  probability  of  the  energy-reducing  (1-2)  collision  occuring 
within  the  time  interval  dt  =  dR/vji(R )  about  R(t)  of  (1-3)  trajectory,  so  that 


dR 

R 


where  Ri  is  the  orbit’s  pericenter  at  which  R  =  0.  On  reversing  the  order  of 
integration, 

OO 

t tcc{E )  —  2tt  J  W(R,  E)  dR 

Ri(0)  o 

where  p\(R,E)  =  jR2[1  —  V(R)/E]  gives  the  maximum  impact  parameter  p\  that 
can  just  acces  the  point  R.  For  purely  atractive  interactions  Ri(0)  =  0  and  V  =  E 
at  Ri(0)  when  repulsion  occurs  at  zero  impact  parameter.  Since 

J  d?  _  2R2  Vu{R) 
l  R(R,E)  v(R)[  E 


then 


oCc{E )  =  47T  J 
Rii  o) 


W(R,E). 

v(R)  1 


E 


dR 


The  (1-2)  collision  probability  within  an  element  ds  of  the  (1-3)  orbit  is 


W(R,  E)dt=  —  =  Nea12(R)v{R)  dt 

A 


where  A  is  the  free  path  length  towards  collision.  Hence  (16)  reduces  exactly  to 


Ro 


(?cc{E)  =  47r  Ne 


VM\ 

E  ) 


R2dR 


which  with  (15)  and  V13  (R)  —  —Ze2/R,  yields 


aCC(E  =  ekT )  =  -ir2ZzRzeNe 


1  +  3(eo  +  3/2)/ e 
2(cq  +  3/2)3 


(17) 


The  factor  ZzRze  originates  from  the  reaction  volume  R%  while  the  additional  R2e 
originates  from  the  (1-2)  scattering  cross  section  (15). 

The  associated  collisional  capture  rate 


OO 

acc  =  (ve)  J  aCc(e)ee_e  de 

0 


is  then 

accM  =  ^Z\v,)I?eN, 

The  rate  for  production  of  all  levels  <  —kTe  is  therefore 

aoc(l)  =  2.2910  -(™)“at.  (18) 

which  is  in  excellent  agreement  with  the  Monte-Carlo  result  of  Mansbach  and  Keck 
[5]  and  a  factor  of  3  lower  than  the  detailed  balance  rate  (12)  derived  from  the 
Thompson  ionization  cross  section  (8).  This  simple  result  ~  ZzT~ 4  5  is  preserved 
even  in  elaborate  classical  treatments  of  the  ionization  [6]  or  of  the  recombination 
[7]  with  correct  averaging  of  the  dynamical  quantities  . 

The  frequencies  vCc  =  otCcNe  and  vRR  =  aRRNe  are  compared  in  Table  (1)  at 
various  Te  for  typical  Ne  ~  108  cm-3.  Since  vcc/^rr  ~  Ne/T{ it  is  seen  that 
Collisional  Capture  dominates  Radiative  Recombination  at  cryogenic  and  lower  Te 
by  large  factors  ~  107.  Also  CC  ~  RR  at  Te  ~  300K  while  CC  <C  RR  at  leV. 


3  e0  + 11/6 
2  (cq  +  3/2)3 


TABLE  1.  Collisional  Capture  and  Radiative  Recombina¬ 
tion  frequencies  for  different  temperatures 


Te 

Vcc{s  J) 

^fift(s_1) 

4  K 

6.3xl04 

5.8xl0-3 

300  K 

2.3X10"4 

3.8x10-“ 

1  eV 

2.0xl0-u 

3.9  xlO"5 

RADIATIVE  DECAY  RATES 


Level  i  decays  radiatively  to  level  j  at  a  rate 

A(i  ->  j)  =  i ME<j/I«)*(S»/9i) 

where  the  line  strength  Sy  in  a.u.  is  J2a,a'  \  (nlm\r\n'l'm')\2  and  A0  =  a3/ 10  = 
1.604  1010  s-1  where  t0  is  the  a.u.  of  time  (e2/a0)/h.  Level  nl  decays  to  all  lower 
levels  s  at  frequency 

ur(nl )  =  ^2  —>  s,l  +  1) 

s<n 

_2  4  1  /  n-5  ,for  1=  n- 1 

3  0  [n®(I  +  \)2  \  4/9  n~ 3  ,  for  l  =  1 

The  circular  (/  =  n  —  1)  states,  mainly  populated  by  collisional  capture  (CC), 
are  therefore  very  long  lived  and  decay  at  frequency 

1010  4  1010 

vr{n,  l  =  n  -  1)  »  — r  <  vr(np  Is)  ^  (s_1) 

Ti  y  Tv 

very  much  less  than  that  for  the  (np  — *  Is)  transition.  The  levels  n0  ~  200  popu¬ 
lated  by  CC  at  4K  have  therefore  a  long  radiative  lifetime  rr  ~  32  s  when  compared 
with  2  10-3  for  the  (n0p  — ►  Is)  transition.  The  overall  rate  of  recombination  via 
capture  into  level  n  followed  by  radiative  decay  is 

vcvT 

ct  ~ - 

Vc  +  Vr 

where  vq  =  vcc  +  vrr-  Since  the  radiative  decay  frequencies  for  n  «  no  are  much 
less  than  the  frequency  accN2  for  collisional  capture  at  ultralow  temperatures, 
the  overall  recombination  rate  is  controlled  by  the  limiting  rate  vr.  The  question 
now  addresed  is  the  further  effect  of  collisions  between  Rydberg  atoms  and  ambient 
electrons  and  ions. 


INELASTIC  COLLISIONS  AND  STARK  MIXING 


Three  important  frequencies  or  timescales  in  the  collision,  depicted  in  Figure  (1) 
are: 

(A)  The  projectile  Rotation  (Collision)  Frequency  ojr  of  the  Projectile 


•  _  bv  large  b  v 
$=  F?  '  b 


1 

Tcoll 


(B)  The  Transition  (Orbital)  Frequency  u>n  of  the  Rydberg  electron: 


FIGURE  1.  The  three  important  frequencies  in  the  Rydberg  projectile  collision 


which  for  transitions  n  ->  n  -  1  between  neightboring  levels  is  simply  u;n,n_i  = 
uo/n3  =  u>n  =  vn/an,  the  orbital  frequency,  where  a„  =  n2a0/Z  and  vn  =  Zv0/n. 

(C)  The  Stark  Precession  Frequency  u >s  for  precession  of  A  about  R: 
Under  the  ion  Electric  Field  S  =  Z\e/R2  the  Stark  frequency 


3  /  c\  3 

Us  ~  2  M(C£)  ~  2 


Z\  CLnVn 
R2 


provides  the  precessional  frequency  of  the  Runge-Lenz  (eccentricity)  vector  A  of 
the  Rydberg  orbit  about  the  field  direction. 

By  considering  the  exp(iu;t)  factor  in  time  dependent  perturbation  theory,  the 
following  types  of  collisions  can  be  characterized  as  in  [8]. 

1.  Orbital  Adiabatic:  ujr  u>n 


CLn  Vn 


where  the  orbital  electron  adjust  itself  adiabatically  to  the  ion  perturbation  and  no 
transitions  occur. 

2.  Orbital  Sudden:  u>r  »  un 

b  v 

—  «  — 
an 


which  causes  inelastic  (n  —*  n ')  transitions 
3.  Stark  Adiabatic:  ujr  <C  u>s 


< 


3  Z\anvn 


bv 

R2 


2  R? 


6  3  ZlVn 

—  «  5  -1— 

where  elastic  collisions  occur. 

4-  Stark  Sudden:  uR  »  us 

b  3  Zivn 

—  >  - - 

an  2  v 

which  causes  /-mixing  ( nl  — >  nl')  transitions. 

5.  The  weak  Field  Condition  us  <  un  i.e.,  the  Stark  splitting  hu>s  <  A 
implies  that  e£  <C  2/3 (e2/a2)  so  that  the  collision  occurs  for  distant  encounters 
at  R  »  an,  the  Rydberg  radius,  for  fields  e£  =  Z\e2 /R2.  The  electron’s  orbital 
time  is  then  much  shorter  than  any  characteristic  time  to  cause  changes  to  the 
elliptical  orbit.  The  vectors  A  and  angular  momentum  L  which  are  constant  for 
the  unperturbed  motion  become  good  dynamical  variables  for  the  description  of 
the  perturbed  motion. 

With  respect  to  orbital  motion,  the  collision  is  sudden  or  adiabatic  according 
to  b  <  be  and  b  >  be,  respectively,  where  be  =  {v/vn)an.  With  respect  to  the 
Stark  frequency,  collision  is  adiabatic  or  sudden  according  to  6  <  6s  and  b  >  bs, 
respectively,  where  bs  =  (vn/v)an.  The  impact  parameter  6-space  can  then  be 
partitioned  according  to: 


Orbital  Sudden 


6C  —  „vn  an 


Orbital  Adiabatic 


Stark  Adiabatic 


(Xr 


Stark  Sudden 


u  3  Z,vn 

us  —  2  t>  Un 


As  v  decreases,  bs  increases  outward  and  be  increases  inward,  thereby  limiting 
the  extent  of  the  sudden  regions  where  n  changes  and  /-mixing  occurs.  The  (v,  6)- 
phase  space  can  be  partitioned  into  the  four  characteristic  regions  illustrated  in 
Figure  (2).  For  v  >  vn,  the  (n, /)  changing  and  /  changing  (Orbital  and  Stark 
Sudden)  shaded  regions  overlap  and  expand,  in  direct  contrast  to  ultracold  speeds 
v  <  vn  where  the  Orbital  and  Stark  Adiabatic  (clear)  regions  increase  and  the 
shaded  regions  diminish  and  do  not  overlap  indicating  few  collisional  changes. 


CLASSICAL  THEORY  OF  STARK  MIXING 

COLLISIONS 


In  addition  to  the  energy  En  and  angular  momentum  L  of  an  unperturbed 
Rydberg  electron  moving  with  velocity  v  in  an  elliptical  orbit  with  eccentricity 


>f$j2  b/a» 


FIGURE  2.  Partitioning  the  v-b  phase  space  map  into  4  regions  characterized  mainly  by: 

(a)  energy  changes,  (b)  energy  and  angular  momentum  changes,  (c)  angular  momentum  changes 
and  (d)  no  changes  under  the  weak  field  condition  (5) 


e  =  (1  —  L2/mank)1^2  and  semi-major  axis  an  =  k/2\En\  in  the  field  —k/r,  the 
Runge-Lenz  (or  eccentricity)  vector 


A  =  d  '[vxL  -  -r] 
r 

vn  =  {v2)x/2  =  k  =  Ze2  u0  =  4~ 

n  n 

directed  toward  the  pericenter  and  normalized  to  angular  momentum  units  is  also 
conserved  in  magnitude  A  —  (mank)l^2e  and  direction.  Moreover  A  •  L  =  0  and 
A2  +  L2  =  mank  =  n2h2.  The  dipole  classically  averaged  over  the  orbit  time  is 


d  = 


3e  A 
2  mvn 


In  a.u.,  A2  =  n2  —  (l  +  1/2)2  and  the  dipole  d  =  (3/2)nA  is  maximum  3/2 n2  for 
highly  eccentric  orbits  l  ~  0  and  is  vanishingly  small  for  circular  orbits  (/  ^  n). 
In  presence  of  an  electric  field  of  intensity  £,  the  angular  momentum  L  changes  at 
the  rate 


dL 

dt 


—er  x  £ 


so  that  the  change  of  L  (classically)  averaged  over  one  orbital  period  T  is  given  by 


Approximation  1:  Assume  that  the  perturbation  is  adiabatic  (u>r  <C  u)n)  with 
respect  to  the  orbital  frequency  a in.  Then  (r  x  £)  =  (r)  x  (£)  so  that 


AL 

T 


— oo  s  x  A 


in  terms  of  the  Stark  frequency 


Ws(t)  =  | ^ 
2  mv„ 


which  is  3n{£au)/2  in  a.u. 

Approximation  2:  Assume  also  that  (£  x  L)  =  (£ )  x  (L)  for  slowly  varying  £. 
It  can  be  shown  that  the  change  A  over  one  period  is 


AA 
•  T 


-<$- 


-u>s  x  L 


Approximation  3:  Replace  the  changes  AL  and  AA  over  one  period  by  LT  and 
AT  to  give  the  set  of  coupled  equations 


dA 

dt 


=  -u>s(t)  X  L 


dL 

dt 


=  -u>s(t)  x  A 


where  u>s  tracks  the  time  variation  of  the  electric  field  in  magnitude  and  direction. 
Under  the  substitution 


L+A  L - A 

2  2 


(19) 


the  above  set  of  differential  equations  becomes  decoupled  to  yield 


=  -vs{t)  x  X  =  +us{t)  x  Y  (20) 

at  flc 

where  X2  =  Y2  =  (L2  +  A2)/ 4  =  n2h2/A.  The  classical  analysis  for  constant 
electric  fields  is  given  Born  in  [9].  For  time  independent  u>s>  both  X  and  Y  precess 
with  constant  frequency  us  about  the  (fixed)  direction  of  u>s-  For  general  time- 
varying  u>s,  the  system  of  differential  equations  (20)  does  not  have  an  exact  solution. 
Percival  and  Richards  [8]  have  used  classical  perturbation  theory  to  solve  (20)  and 
then  provided  a  diffusional  theory  of  angular  momentum  mixing.  Bellomo  et  al 
[10]  approached  the  same  problem  by  proceeding  via  the  time  evolution  propagator 
U±ot(t,  t0)  for  X  and  Y  in  the  rotating  frame,  an  approach  which  results  in  formulae 
too  complicated  for  physical  changes  AL  and  AA  to  be  extracted. 

An  Exact  Analytical  Solution  is  however  possible  when  the  external  time 
dependence  of  <us{t)  is  provided  by  the  electric  field  £  =  Zxe(—R)/ R2  of  a  projectile 
ion  or  electron  of  charge  Zxe  passing  the  Rydberg  atom  at  large  distances  R  »  an, 
the  condition  for  the  weak  field  approximation.  Then 


which  varies  in  both  magnitude  and  direction.  Since  the  vector  R  rotates  at  the 
frequency 

d$  „  vb  „ 

Ws(()  =  - Tt 1  =  W) 1 


then  the  vector 


/.X  ws 

a{t)  =  — 

Ur 


R(t) 

2  \Z  )  bv  w 


has  the  constant  magnitude  3/2(Zl/Z)(bv)~1  for  a  particle  with  reduced  speed  v  = 
v/vn  in  a  given  trajectory  of  reduced  impact  parameter  b  =  b/an.  Since  the  value 
a  =  l  separates  the  Stark  Sudden  (a  <C  1)  and  the  Stark  Adiabatic  (a  2>  1)  regions 
at  b  =  bs  (see  Figure  (2)),  we  refer  to  a  as  the  Stark  parameter.  This  suggests 
that  $  could  be  a  more  useful  variable  instead  of  time  t.  Hence  the  first  equation 
in  (20)  is 


dX 


dXd$ 
dt  dt 


a(0,  sin^cos'F)  x  X 


(21) 


Since  the  tip  of  X  moves  on  the  surface  of  a  sphere,  the  solution  of  this  equation, 
the  final  X  vector  (corresponding  to  t  — ►  oo  or  <f>  =  0),  is  a  rotation  of  the  initial 


X  vector  (corresponding  to  t  — *■  — oo  or  $  =  ir).  The  second  equation  in  (20)  has 
a  similar  solution.  The  exact  form  of  the  resulting  rotation  matrices  are  presented 
in  the  Appendix.  As  a  result  of  the  collision  the  vectors  L  and  A  change  to 

L\  i  1  +  a2cos(7T7)Li  — a27sin(7T7)A2  a(— 1  +  cos(7T7))A3 

L2  —  ~  72cos(tt7)L2  7sin(7r7)L3  — a7sin(7T7)Ai  (22) 

L3  ^  7sin(7r7)L2  — (a2  +  cos(7T7))L3  a(l  —  cos(7T7'))Ai 

and 

Ax'  ^  ’  (1  +  a2  cos(7T7))Ai  — a27sin(7T7)L2  a(-l  + cos(7T7))L3 

A2  =—  — 72  cos(7T7)A2  7sin(7T7)A3  a7sin(7T7)Z,i 

A3  ^  7sin(7T7)A2  —  (a2  +  cos(7T7))A3  a(l  —  cos(7r7))L! 

(23) 

where  7  =  \/l  +  a2  and  the  components  of  the  initial  and  final  vectors  (L,  A)  and 
(L,  A)  are  defined  in  the  fixed  coordinate  frame  of  Figure  (1).  It  easy  to  check 
that  the  above  exact  solutions  satisfy  the  invariant  relations 

L  •  A  =  L  •  A  =  0 

and 

L2  +  A2  =  L2  +  A2  =  n2h2 

The  new  (n,  L)  orbit  is  confined  to  a  plane  perpendicular  to  the  new  L  and  the 
energy  is  preserved  (n  is  not  changed).  The  key  to  the  present  exact  solutions  of 
(20)  for  a  non-uniform  rotating  field  is  provided  by  the  recognition  that  $  is  a  more 
useful  variable  than  time  t  and  by  the  Rotation  Matrix  Algebra  (in  the  Appendix). 

In  the  Stark  Sudden  region,  the  Stark  parameter  o;  is  <C  1.  In  this  limit,  the 
solutions  have  the  simpler  form 

Lx  (1  —  2a2)  Lx  —  2aA3 

~L2  =  L2  +  it  a2/ 2  1 L3 

Z3  _  (1  -  2 a2)L3  -  7ra2/2  L2  +  2aAx  _ 

and 

'Ax  1  [  (1  -  2a2)Ai  -  2aLz 

A2  =  A2  +  ira2 /2  A3 

_  A3  _  (1  -  2a2) Az  -  7ra2/2  A2  +  2aLx  _ 

The  magnitude  of  the  final  angular  momentum  (for  a  <C  1)  is  then 

L2  =  L2  +  4a(— A$Lx  +  AiL3)  +  4a2(A2  —  A2  —  L2  +  L2) 

On  retaining  only  the  first  order  term  in  a 

L2  «  L2  +  4a (A  x  L)  •  y  =  L2  +  4aALcos^ 


where  £  is  the  angle  between  A  x  L  and  the  Oy  axis  (Figure  (1)). 

Since  the  vectors  A  and  L  are  always  orthogonal,  the  angle  £  is  then  uniformly 
distributed  over  the  range  (0,7r)  if  L  is  uniformly  distributed  on  the  sphere  with 
radius  L,  as  for  an  equilibrium  distribution  of  magnetic  m-sublevels.  The  initial 
L-states  then  evolves  to  the  final  L  states  with  angular  momentum  distributed 
according  to 


T(L)  dL  =  —  d£  = 

7T 


2 _ L _ 

*  y/(4aAL)2  -  (L2  -  L2)2 


dL 


In  fact  T  can  be  identified  with  WL_i,  the  transition  probability  density  for  L  — *•  L 
transitions. 

For  arbitrary  a,  the  transition  probability  W  has  in  general  not  a  simple  form.  It 
is  however  important  to  recognize  that  the  probability  W  depends  on  the  projectile 
velocity  v  and  impact  parameter  b  only  via  the  Stark  parameter  a  =  2>Z\n/2Zvb. 

The  cross  section  for  /-mixing  in  the  the  Rydberg  atom  is  now 


da 

dL 


1= 2-7  Mb  - 2-  (I)2  (§)2  /  dQ 

bs  0 


where  the  a-integral  is  independent  of  the  speed  of  the  projectile  and  the  Rydberg 
level  n.  This  general  result  (where  I  is  the  above  mentioned  integral) 


is  in  agreement  with  the  behavior  of  a  normalized  Born  approximation  presented 
in  the  recent  review  [11]. 


SUMMARY 

In  addition  to  providing  radiative  and  collisional  rates,  the  three  main  points  of 
this  paper  are: 

•  Development  of  a  simple  classical  expression  for  the  cross  section  (16)  and  rate 
(18)  of  the  three-body  collisional  capture  (CC)  into  level  n.  The  result  agrees 
with  the  Monte-Carlo  trajectory  calculations  [5]  and  preserves  the  Z3T~ 4  5 
dependence  of  a. 

•  Presentation  of  an  exact  classical  analytical  solution  for  the  changes  AL  and 
AA  in  the  angular  momentum  L  and  Runge-Lenz  vector  A  of  a  Rydberg  atom 
in  the  time  varying  electric  field  £  =  —Z\eR{t)/R2(t). 


•  Pointing  out  that  recombination  at  ultracold  temperatures  proceeds  by  very 
rapid  three-body  collisional  capture  into  levels  n/,  with  high  /,  followed  by 
angular  momentum  mixing  by  collisions  mainly  with  ions.  Recombination 
then  becomes  stabilized  by  radiative  decay  of  the  lower  /-levels  so  produced. 
This  sequence  is  in  direct  contrast  to  recombination  at  Te  >  1  eV  where 
collisional  and  radiative  capture  are  generally  the  rate  limiting  steps. 


APPENDIX:  ROTATION  MATRIX  ALGEBRA 


For  a  given  direction  n,  rotate  vector  x  by  angle  <f>  about  the  fixed  direction  n  to 
give 

x1  =  cos0x  +  (1  —  cos  4>){n  •  x)n  +  sin  0(n  x  x)  =  R(h,  </>)x 
This  rotation  corresponds  to  the  matrix  operator: 


R(h,  4>) 


c  +  (1  —  c)n\  (1  —  c)niri2  —  sn3  sn2  +  (1  —  cri\n3 

(1  -  c)nin2  +  sn3  c  +  (1  -  c)n22  ~sn i  +  (1  -  c)n2n3 
—sn2  +  (1  -  c)nin3  sni  +  (1  -  c)n2n3  c  +  (1  -  c)n32 


where  c  =  cos 0,  s  =  sin 0  and  h  =  h(ni,n2,n3).  This  matrix  can  also  be  written 
as 

R(n,  4>)  =  exp  4>{n\Ix  -1-  n2Iy  +  n3Iz ) 
where  IX)  Iy,  Iz  are  the  infinitesimal  generators  of  the  0(3)  group,  i.e. 


'  0 

0 
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The  equation  (21)  can  now  be  rewritten  as 


dX 

</$ 


=  a(cos  $/2  +  sin  3>/y)X 


The  transformation  (“  rotation  frame”) 


X'  =  R(x,  $)X 


0-10 
1  0  0 

0  0  0 


eliminates  the  $  dependence  of  the  precessional  frequency  of  the  vector  X'  in  the 
rotating  frame  since 


^  =  (^/T1  +aR(cos$iz  +  sin$4)R-1)X'  =  (!x  +  aIz)X' 
The  general  solution  for  X'  is  then 

X'(4>)  =  exp  ((*  -  %)(!„  +  al,j)  X'(4>0) 


or  in  the  language  of  rotation  matrices, 


X'($)  =  R 


1  +  a 


■  (hi  +  cm3),  ($  -  $o 


a2  X'($0) 


The  exact  solution  for  X(t)  in  the  laboratory  frame  Figure  (1)  is 


X($)  =  R(hu-^)R 


'l  +  a2 


{h\  +  an3),  ($  -  $o)V^l  +  a2  R(h i,  <&0)X($0) 


Initially  X  =  X*  when  the  projectile  approaches  from  infinity  so  that  =  7 r.  The 

collision  is  completed  when  j  =  it  so  that 


Xf  =  R 


/l  +  a; 


:(hi  +  ah3),  -7TV1  +  Q:2  i?(ni,7r)Xi 


Upon  replacing  a  by  —a  the  solution  Y/  of  the  second  equation  in  (20)  is 


Yf  =  R 


+  a2 


(hi  -  ah3),  -7rvl  +  a2  R(hi,  7r)Yj 


These  last  two  equations  and  the  substitutions  (19)  provide  the  relations  (22) 
and  (23)  between  the  initial  and  final  angular  momentum  and  Runge-Lenz  vectors. 
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The  classical  atomic  form  factor 
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The  general  trends  exhibited  in  the  variation  of  the  inelastic  form  factor  in  collisiona!  transitions 
nl  — ►  n'f ,  when  l'  is  changed  and  n,  l  and  n!  are  kept  fixed  are  explained  solely  in  terms  of  classical 
mechanics.  Previous  quantal  results  are  reproduced  from  purely  classical  mechanics  principles. 
Our  conclusions  are  valid  not  only  for  large  quantum  numbers  (which  provide  the  usual  classical 
correspondence)  but  also  for  other  cases,  which,  up  to  now  have  only  been  described  by  quantal  or 
semiclassical  methods.  The  interesting  trends  exhibited  in  the  form  factor  are  directly  reflected  in 
experimental  and  theoretical  treatments  of  collisions  involving  excited  atoms. 
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With  the  advent  of  new  technology  which  facilitates 
the  accurate  measurement  [1]  of  electron  -  excited  atom 
collision  cross  section  there  has  also  been  renewed  interest 
in  the  theory  [2]  of  collisions  involving  Rydberg  atoms. 
Recent  experiment  [1],  in  particular,  has  confirmed  that 
the  cross  section  for  the  quadrupole  23S  — >  3 3D  transi¬ 
tion  in  e  —  He(2sS)  collisions  is  much  higher  then  that 
for  the  pure  dipole  23S  — >  33P  transition  at  low  and 
intermediate  energies,  in  accord  with  the  theoretical  pre¬ 
dictions  of  [3]  (Born  and  Multichannel  Eikonal  approxi¬ 
mations).  Flannery  and  McCann  [4]  have  noted  that  this 
unexpected  behavior  is  only  part  of  a  more  general  sys¬ 
tematic  trend  in  that  (a)  the  23S  — >  nzD  collisional  tran¬ 
sitions  are  predominant  over  all  other  transitions  to  the 
same  value  n-value,  even  for  transitions  to  the  electronic 
continuum  and  (b)  that  there  is  a  unique  value  Vmax  of 
the  final  angular  momentum  /'  that  is  preferentially  pop¬ 
ulated  in  nl  — >  nfV  transitions  (n'  n)  in  collisions 
between  Rydberg  atoms  and  electrons  or  atoms. 

The  origin  of  this  general  behavior  was  traced  [4]  to 
the  variation  with  lf  of  the  quantum  mechanical  inelastic 
form  factor 

ffi( q)  =  (^/ (r) I l^i (r ))  =  (<pf{ p  +  q)|0i(p)>  (1) 

for  i(n,l)  —»  f(nf,V)  transitions  between  atomic  states; 
tpij{ r)  are  the  wave  functions  in  position  space  and 
p)  —  (2?r^)-3/2  f  ipij(r)  exp(-zqr/ft)dr,  the  wave 
functions  in  momentum  space. 

When  an  instantaneous  impulse  applied  at  t  =  to 
transfers  momentum  q  to  an  atomic  electron,  the  exact 
solution  of  Schrodinger’s  equation  under  Hamiltonian 

H( p,  r,  t)  =  p2 /2m  —  e2/r  ~  r  •  q 8(t  —  to)  (2) 
is 

^(r,  t)  =  (1  +  eiqr/h0(t  -  t0))ipnim{r) 

where  6  is  the  Heaviside  step  function.  The  probability 
for  i  =  | nl)  — ►  /  =  |n'/')  transitions  from  the  (21  4-  1) 
initial  sublevels  is  then 


Pnl,n’l’(q)  =  \Wn’l'm\2=  £  |<nTm'|ei^/'V™>|2 

m,m' 

(3) 

The  probability  of  any  impulsive  i  — ►  /  transition, 
whether  due  to  particle  collisions  or  electromagnetic  field, 
is  therefore 

Pu(  q)  =  l^/.(q)l2  (4) 

which  provides  physical  significance  to  the  inelastic  form 
factor,  a  fundamental  property  of  the  atom.  For  impul¬ 
sive  collisions  between  a  particle  1  and  a  Rydberg  elec¬ 
tron  2  bound  to  a  core  3,  the  overall  transition  matrix 
element  T  decomposes  as  [5] 

TiM)  =  J),(q)T12(q)  (5) 

where  Ti2,  the  matrix  element  for  (1-2)  free- free  elastic 
scattering  in  the  (1-2)  center  of  mass,  is  a  function  only 
of  q,  as  for  Coulomb  scattering  T12  =  47t h2e2/q2,  or  for 
Born’s  approximation,  Ti2  =  f  V(ri2)exp(iqr/h)dri2- 
The  probability  of  transition  in  the  target  atom  per  each 
(1-2)  impulsive  encounter  is  P*/  =  |Tj/|2/|Ti2|2,  in  agree¬ 
ment  with  (4). 

The  cross  section  is  obtained  by  the  following  integra¬ 
tion  of  the  form  factor  (5)  over  momentum  change, 

ki+kf 

Oif  =  §  (M/M12)2  J  \rfi(q')\2  l/i2(</)|2  q'dq'  (6) 

\ki-kf\ 

where  kij  are  the  initial  and  final  wavenumbers  of  rela¬ 
tive  motion  of  the  projectile  -  target  system  of  reduced 
mass  M  and  q(  =  h\k{  —  k/|  is  the  momentum  change. 
The  scattering  amplitude  for  (1-2)  collisions  of  reduced 
mass  M\2  is  /12  =  (2Mi2/47rft2)Ti2.  For  (1-2)  slow  col¬ 
lisions  with  scattering  length  a,  the  Fermi  interaction 
y(r12)  =  [47ra(ft2/M12)]6(ri  —  r2)  also  yields  decomposi¬ 
tion  (5)  with  fi2  —  a . 
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'  .  The*  inelastic  quantal  form  factor  therefore  not  only  ex¬ 
erts  primary  importance  in  collision  studies,  but  also  has 
a  deep  physical  reality.  In  recent  experimental  studies 
of  excitation  of  Rydberg  atoms  by  short  unipolar  half¬ 
cycle  electromagnetic  pulses  the  transition  amplitude  is 
determined  directly  by  the  inelastic  form  factor  [6] . 

Analytical  quantal  [7,8]  and  semiclassical  [9]  form  fac¬ 
tors  are  available  although  general  systematic  trends  can 
not  be  easily  extracted  from  them.  A  key  point  of  this  pa¬ 
per  is  that  a  complementary  classical  approach  can  also 
be  developed  in  a  way  which  reveals,  quite  succintly,  im¬ 
portant  aspects  which  remain  hidden  within  the  quantal 
treatment. 

Consider  a  Rydberg  atom  in  a  stationary  (n,  /)  state 
with  energy  E  and  angular  momentum  L.  If  the  atom  is 
perturbed  by  any  general  impulsive  field  (as  in  eq.  (2) 
or  the  Fermi  interaction),  then  the  transition  probabil¬ 
ity  to  the  final  state  (n',  V)  (of  energy  E'  and  angular 
momentum  Lf)  is  the  inelastic  form  factor. 

The  quantal  probability  density  for  finding  the  electron 
in  the  radial  interval  (r,  r  +  dr)  is 

Pqni(r)  =  r2\Rni\2  (7) 

where  Rni  is  the  hydrogenic  radial  wave  function  ex¬ 
pressed  in  terms  of  the  generalized  Laguerre  polynomial. 

The  phase  space  of  a  classical  atom,  with  Hamiltonian 
H( r,  p)  =  p2 /2m  4-  V(r),  angular  momentum  L(r,  p)  = 
rxp  and  period  t„j  in  stationary  state  (n,  l)  is  populated 
according  to  the  microcanonical  distribution  [10] 

Pni  drdp  =  {(2nh2/rnl)d(H  -  £)<5(|L|  -  L)} 

(8) 

normalized  to  (2/  +  1)  states  in  all  of  phase  space.  On 
integrating  over  the  momentum  space  p  and  angular  part 
of  the  configuration  space  r,  the  classical  distribution  is 

Pni(r)  dr  = -  T  dr 

Tnl  r 


R±  =  n2  {1  ±  e}  =  n2  jl  ±  [l  -  (/  +  l/2)2/rz2] 1/2  } 

(10) 

The  classical  distribution  is  zero  outside  the  accessible 
region,  bounded  by  R±. 

By  using  definition  (1)  the  transition  probability  (3) 
can  be  converted  to  the  new  form 

^i/(q)  =  (2nh)3  j Pni(r»P)Pn'i'(r,P  +  q)  drdp  (11) 

where  the  quantal  distributions  in  phase  space  are  given 
by  pq  =  (27rft)"’3/2^(r)exp(-zp  •  r/h)<f>*(p).  This  form 
is  now  suitable  for  classical  correspondence  on  replacing 
densities  pq  by  the  phase  space  distributions  (8).  The 
basic  definition  of  the  classical  form  factor  is  therefore 
given  by  (8)  and  (11).  The  physical  significance  is  that 
the  initial  and  final  state  correspond  to  definite  regions  in 
phase  space,  populated  according  to  the  microcanonical 
distribution  (8),  and  the  transition  probability  is  given, 
in  a  geometric  sense,  by  the  amount  of  overlap  of  these 
regions.  In  configuration  space,  the  regions  are  spherical 
shells  with  inner  and  outer  radii  given  by  eq.  (10),  the 
pericenter  (Jl”)  and  apocenter  (R+)  of  the  Kepler  orbit. 

Analytical  expressions  with  explicit  dependence  on  q 
for  quantal  and  classical  transition  probabilities  are  de¬ 
veloped  in  a  separate  paper  [8].  Rather  than  examining 
the  V  variation  of  (11)  for  a  given  q ,  the  key  results  are 
more  readily  deduced  and  are  easily  transparent  by  in¬ 
vestigating  the  probability  for  all  momentum  transfers 

Fnl-*n'V  =  J  pi/(q)  dq  =  (27Tfi)3  J  dr  pcnl(r)  pcn,  v  (r) 

TZ 

(12) 

where  7Z  is  the  overlapping  region  in  configuration  space 
defined  by  intersection  of  (R~,R+)  and  (R7,Rt)  inter¬ 
vals. 

Inserting  p(r)  =  4npc(r)r2  with  (9)  in  (12)  gives  the 
classical  form  factor  (CFF) 


where  the  radial  speed  is  given  by  rar2/ 2  =  E-V(r)  — 
(Z-bl/2)27i2/2mr2.  For  the  Kepler  atom  (rn/  =  27rn3  a.u) 
and  pcnl  (in  a.u.)  is 

2 _ 1 _ ( l  4- 1/2)21 -1/2  _  J_  J_ 

r  n2  r2  ttti^  r(r) 

(9) 

The  quantal  (7)  and  classical  (9)  radial  probability 
densities  are  illustrated  in  figure  Fig.  1.  As  in  the  text¬ 
book  example  of  the  harmonic  oscillator,  the  classical 
distribution  has  singularities  at  the  corresponding  turn¬ 
ing  points  given  by  the  radii  (in  a.u.): 


PCn,(r)  =  -T 


7T  71° 


fim  a  x 

(2/  + 1)  7  dr/r^ 

n3n/3  J 

f^min 


(13) 


where  Rmin  —  max(R^  ,Ry)  and  Rmax  —  min(R^,J2j) 
define  the  bounds  of  the  overlapping  region  7 Z.  Different 
overlap  situations  are  illustrated  in  Fig.  2  for  a  repre¬ 
sentative  case.  The  gray  region  is  the  accessible  region 
for  the  initial  state  and  the  curves  are  possible  final  state 
trajectories.  Transitions  only  occur  when  the  final  state 
trajectory  penetrates  the  initial  state  accessible  region. 
The  longer  time  spent  by  the  electron  on  the  final  state 
trajectory  within  the  initial  state  accessible  region,  the 
bigger  is  the  transition  probability. 
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.  As  l  increases  from  zero  to  its  maximum  value  for  cir¬ 
cular  orbits,  R~  increases  from  zero  to  n2,  while  R+  de¬ 
creases  from  2n2  to  the  same  value  n2.  For  final  states 
n'  >  \f2n,  then  Rmax  =  Rt  f°r  values  of  Z'.  Three 
regions  of  overlap  are  then  apparent  and  are  respectively 
accessed  as  Z'  is  increased: 

•  Region  I,  RJ  <  RJ  :  Here  the  overlap  region 
1 Z  ~  (Rj,Rf)  is  determined  solely  by  the  initial 
state  and  has  spatial  extent  which  remains  constant 
as  Z'  is  varied  from  zero  to  some  value  Zi  where 
RJ  =  RJ .  There  is  always  an  orientation  of  the 
final  orbit  which  will  then  intersect  the  initial  or¬ 
bit,  as  exhibited  in  Fig.  2,  for  (n  =  3,  1  =  2)  and 
(n'  =  8,  V  =  0  —  2)  orbits.  The  l '  variation  of  (13) 
is  contained  solely  within  the  increasing  integrand 

•  Region  II,  RJ  <  RJ  :  Here  the  overlap  region  7 Z  = 
(RJ,  Rj)  includes  the  f-pericenter  and  has  spatial 
extent  which  decreases,  as  V  increases,  eventually 
to  zero  when  RJ  =  Rf.  In  this  region,  the  initial 
and  final  orbits  can  intersect  each  other,  as  for  the 
(n'  =  8,  V  =  4)  orbit  in  Fig.  2.  The  V  variation 
of  (13)  results  from  variation  of  both  the  increasing 
lower  limit  RJ  and  the  increasing  integrand  (r/)”1. 

•  Region  III,  RJ  >  Rf  :  Here  the  initial  and  final 
trajectories  no  longer  intersect,  since  the  pericenter 
of  the  final  state  is  greater  that  the  apocenter  of  the 
initial  state.  This  region  where  (i  — >  /)  transitions 
do  not  occur,  as  illustrated  by  (n'  =  8,  Z'  =  5, 6, 7) 
orbits  in  Fig.  2,  is  the  classically  inaccessible  re¬ 
gion. 

The  boundaries  between  regions  I  and  II  and  between 
regions  II  and  III  occur  respectively  at  V  =  l\  where 
RJ(nr,V)  =  RJ(n,l)  and  at  Z'  =  I2  where  RJ(n',l')  = 
Rf(n,l).  Thus  l\  and  Z2  are  given  by 

(h,2  +  =  n2(  1  T  e)  (2  —  (1  =f  e>2/n'2)  (14) 

where  e  is  the  eccentricity  [1  —  (Z  +  1/2 )2/n2]1/2  of  the 
initial  orbit. 

Variation  of  the  CFF  (13),  with  final  angular  momen¬ 
tum  Z '  is  then  determined  both  by  the  lower  integration 
limit  Rmm  (which  is  a  constant  RJ  in  region  I  and  in¬ 
creases  as  RJ  in  region  II)  and  by  the  integrand  (r/)_1. 
Fig.  3  illustrates  the  general  pattern.  As  l '  is  increased 
from  0  to  Zi  (region  I),  increase  in  CFF  originates  purely 
from  the  increasing  integrand  (r/)”1.  As  l'  is  varied 
from  li  to  I2 ,  the  increasing  integrand  is  offset  by  the 
decreasing  range  (RJ ,  Rf)  of  integration  (region  II).  For 
I2  <  V  <  n  —  1,  CFF  is  zero  because  transitions  are  not 
classically  allowed  in  region  III. 


At  /'  =  Zi  the  trajectories  touch  only  at  their  corre¬ 
sponding  pericenters  and  CFF  has  a  turning  point  sin¬ 
gularity  characteristic  of  classical  descriptions.  The  zero 
radial  speed  of  the  electron  at  the  contact  point  of  both 
initial  and  final  orbits  causes  the  infinite  CFF  (transition 
probability). 

As  is  evident  from  Figs.  3-5,  the  agreement  between 
the  classical  and  quantal  results  is  excellent  in  region  I, 
even  for  small  quantum  numbers.  In  region  II,  the  quan¬ 
tal  results  oscillate  about  CFF.  Since  classical  motion  is 
confined  to  a  definite  region,  the  dramatic  fall  for  large  l' 
is  more  steep  than  that  for  the  quantal  case  where  states 
have  exponential  tails  within  the  classical  inaccessible  re¬ 
gion  III.  As  expected  from  correspondence  principles,  for 
the  larger  quantum  numbers,  the  quantum  form  factor 
tends  to  CFF,  even  in  the  regions  II  and  III,  as  shown 
in  Fig.  4.  The  quantal  results  exhibit  maxima  in  the 
neighborhood  of  /'  =  Zi ,  Z2  where  CFF  has  the  classical 
singularities.  The  position  of  l\  defined  by  (14)  in  the 
limit  of  large  Z,  where  the  eccentricity  e  — >  0,  is 

Zi(Z  — >  n  -  1)  =  nV 2(1  —  1/2  (n/n')2)1/2  -  1/2 

an  exquisite  result  for  initial  circular  orbits.  For  n'  »  n, 

1 1  tends  from  the  bottom  to 

Zi(Z  — ►  n  -  l,n'  >  n)  =  n\Pl  -  1/2  (15) 

a  key  result  in  detailed  agreement  with  that  previously 
derived  from  consideration  of  the  quantal  momentum  - 
space  overlap  [4]. 

For  small  initial  angular  momentum  1,  e  — ►  1  and  l\  is 
then  zero  so  that  the  maximum  CFF  is  given  by 

l2(l  —  0)  =  2n(l  -  (n/n')2)1/2  -  1/2 

appropriate  to  highly  eccentric  initial  orbits.  In  the  n' 
n  limit  then 

Z2(Z  — *  0,n'  n)  =  2n  —  1/2  (16) 

As  the  initial  1  increases,  there  is  therefore  a  slow  varia¬ 
tion  (2n  — *  >/2n)  in  the  position  of  the  maximum  of  CFF 
(Eq.  (15)  and  (16)),  which  is  pushed  slightly  to  lower  val¬ 
ues.  This  theoretical  prediction  is  also  confirmed  by  the 
quantal  results  [4]. 

When  the  energy  Ef  of  the  final  orbit  is  not  sufficient  to 
accommodate  the  value  of  Z2  deduced  above  (nr  <  \/2n), 
the  peak  in  CFF  (as  in  Fig.  5)  is  given  by  Zi,  provided  the 
initial  Z  is  large  enough.  When  the  final  n'  is  sufficiently 
small  so  that  the  lower  l\  can  not  be  accomodated,  i.e. 
1 1  >  nf  —  1  CFF  and  the  quantal  result  exhibit  a  mono¬ 
tonic  increase  confined  to  Region  I,  which  is  always  char¬ 
acterized  by  excellent  agreement  between  quantal  and 
classical  results.  Classical  expressions  (11,13)  for  the 
atomic  form  factor  have  been  derived.  The  pattern  ex¬ 
hibited  by  the  Z'  variations  (Figs.  3-5)  is  essentially  iden¬ 
tical  with  the  quantal  pattern.  The  positions  of  maxima 
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'  of  the  ? '-variation  of  CFF  depend  strongly  on  the  initial 
n  and  only  weakly  on  the  initial  Z,  in  agreement  with  the 
quantal  calculations  [4],  which  were  restricted  to  certain 
cases.  Excellent  quantitative  agreement  between  classi¬ 
cal  and  quantal  results  makes  the  classical  form  factor  a 
very  useful  tool  particularly  at  large  quantum  numbers 
(Rydberg  atoms)  where  exact  quantal  results  are  not  easy 
to  obtain  (either  analytically  or  numerically)  and  to  use, 
due  to  the  highly  oscillatory  nature  of  the  wave  function. 
Although  the  emphasis  here  is  on  the  electron  form  fac¬ 
tors,  the  present  analysis  is  applicable  also  to  form  factors 
for  transitions  between  rovibrational  states  of  molecules. 

This  work  was  supported  by  AFOSR:  F  49620-96-1- 
0142  and  NSF. 


FIG.  2.  Various  final  state  ( n  =8,  l'  —  1  —  7)  tra¬ 
jectories  and  the  initial  accessible  region  corresponding  to 
(n  =  3,  l  =  2). 
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FIG.  3.  Characteristic  dependence  of  the  the  inelastic  form 
factor  on  the  final  angular  momentum  /',  for  fixed  n  (=  3), 
l  (=  2)  and  n'  (=  8).  Classical  calculations:  solid  line;  quantal 
results:  dots. 
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FIG.  4.  Classical  (solid  line)  and  quantal  (dots)  inelastic 
form  factor  for  transitions  from  state  (n  =  35,  l  =  30)  to 
( n'  =  55,  /'  =  0  — ♦  54)  states. 


FIG.  1.  Classical  and  quantal  radial  densities  of  probability 
of  localization  for  the  stationary  state  of  the  hydrogen  atom 
(E=  -1/(2  x  202)  and  Z  =  8). 
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FIG.  5.  Classical  (solid  line)  and  quantal  (dots)  inelastic 
form  factor  for  transitions  from  n  =  6,  /  =  5  to  n'  =  8  states. 
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Cross  sections  for  electron  excitation  out  of  the  2  3S  metastable  level  of  He  into  the  2  P.  3  5.  3  P . 
3  }D.  4  ?S.  4  3P.  4  3D,  5  3S.  and  5  3D  levels  have  been  obtained  for  energies  up  to  18  eV  We  have  observed 
a  broad  excitation  function  for  the  2  :'P  level  with  a  peak  apparent  cross  section  of  1.2x  10' 14  cm*.  For  the 
n  =  3,  4.  and  5  levels,  the  excitation  functions  show  a  pattern  of  sharp  peaks  for  excitation  into  the  n  "S  le\els. 
slightly  less  sharp  peaks  for  excitation  into  the  n  }P  levels,  and  relatively  broad  peaks  for  excitation  into  the 
n  3D  levels.  Absolute  cross  sections  have  been  obtained  for  all  the  above  mentioned  levels  using  a  laser- 
induced  fluorescence  technique  and  the  results  agree  well  with  experimental  values  reported  by  Lagus  et  ai. 
[Phys.  Rev.  A  53,  1505  (1996)]  The  cross  sections  for  the  2'^S^n  :'P  excitations  which  correspond  to 
dipole-allowed  optical  transitions  are  smaller  than  the  corresponding  2  3S— }S  and  2  'S  —  n  D  excitation 
cross  sections,  in  contrast  to  the  trends  observed  for  excitations  out  of  the  ground  level.  This  reversal  behavior 
is  discussed  in  terms  of  the  dipole  matrix  element  sum  rule.  Our  cross-section  data  are  compared  with  those  of 
the  alkali-metal  atoms.  [S  1050-2947(97)06204-5] 

PACS  number(s):  34.80.Dp 


I.  INTRODUCTION 

Excitation  out  of  the  metastable  levels  of  rare  gases  is  an 
important  mechanism  in  a  wide  variety  of  phenomena,  such 
as  high  density  gas  discharges,  astrophysical  plasmas,  and 
electron-beam  pumped  lasers.  In  the  case  of  the  2  lS  and  2  3S 
metastable  levels  of  helium,  much  theoretical  work  has  been 
done  on  calculating  both  differential  and  integral  cross  sec¬ 
tions  for  electron  excitation  out  of  these  excited  levels  [1-7]. 
An  early  measurement  by  Gostev  et  al.  [8]  shows  results  that 
are  in  serious  disagreement  with  the  subsequent  experiment 
of  Mityureva  and  Penkin  [9].  Measuring  cross  sections  for 
processes  out  of  the  metastable  levels  has  proven  to  be  much 
more  difficult  than  measuring  cross  sections  for  the  corre¬ 
sponding  processes  out  of  the  ground  levels  of  these  gases. 
This  is  primarily  due  to  the  difficulty  in  producing  sufficient 
densities  of  metastable  atoms  to  study  and  separating  the 
associated  signal  from  signals  due  to  ground  level  atoms.  In 
our  laboratory  we  have  used  a  hollow. cathode  discharge  to 
produce  a  thermal  atomic  beam  containing  metastable  he¬ 
lium  with  a  density  of  6X107  cm-3.  With  this  apparatus, 
along  with  the  use  of  special  techniques  for  detecting  very 
low  photon  emission  signals,  we  have  measured  excitation 
functions  and  cross  sections  for  electron  excitation  out  of 
both  the  singlet  and  triplet  metastable  levels  for  electron  en¬ 
ergies  up  to  the  onset  of  ground  level  excitation  [10-12]. 
Our  current  work  represents  a  refinement  and  extension  into 
the  n  =2  and  n  -5  triplet  levels  of  the  previous  results  of  Rail 
et  ai  [10].  We  have  found  distinct  patterns  in  the  shapes  of 
the  excitation  functions  for  the  5,  P,  and  D  levels. 

Part  of  the  motivation  for  our  effort  is  that  the  experiment 
of  Rail  et  ai ,  which  used  the  aforementioned  thermal  atomic 
beam  as  a  target,  and  the  experiment  of  Lagus  et  ai  [13], 
which  used  a  fast  (1  keV)  beam  of  metastable  atoms  pro¬ 
duced  by  charge  exchange  collisions  as  a  target,  had  a  dif- 
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ference  in  their  absolute  cross  sections  of  about  5(Fr.  We 
feel  that  it  is  important  to  improve  the  absolute  calibration  of 
the  thermal  atomic  beam  experiment  to  determine  the  source 
of  this  discrepancy.  As  discussed  in  Sec.  Ill  E.  with  the  im¬ 
proved  absolute  calibration  our  absolute  cross-section  mea¬ 
surements  using  the  thermal  metastable  beam  target  show 
good  agreement  with  the  results  of  Lagus  et  al. ,  which  are 
based  on  an  entirely  different  method  of  absolute  calibration. 

II.  EXPERIMENTAL  APPARATUS 

A  schematic  diagram  of  our  experimental  apparatus  is 
shown  in  Fig.  1,  which  is  quite  similar  to  the  apparatus  de¬ 
scribed  by  Lockwood  et  ai  [11.12].  In  this  experiment,  an 
electron-beam  crosses  at  right  angles  an  atomic  beam  con¬ 
taining  both  ground-state  and  metastable  He  atoms  produced 
by  a  hollow  cathode  discharge.  The  emission  intensity  from 


FIG.  1.  Schematic  diagram  of  the  experimental  apparatus  used. 
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